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Abstract 

The paper introduces a semi-analytical method for calculating the Hall conductivity in the single- 
band approximations. The method goes beyond the linear response theory and, thus, imposes no 
limitation on the electric fields magnitude. It is shown that the Hall current decreases with increase 
of the electric field, if the Bloch frequency (which is proportional to the electric field) exceeds the 
cyclotron frequency (which is proportional to the magnetic field). The obtained results can be 
directly applied to the system of cold Bose atoms in a 2D optical lattice, where the single-band 
approximation is well justified. 
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I. INTRODUCTION 

Begging from works by Ohm in XIX century [l| and till seventeenth of XX century all 
studies of ordinary conductivity in solid crystals were restricted by the weak field regime, 
where the electric field is considered as a perturbation. This approach was actually justified 
because for typical laboratory conditions the Bloch frequency, which is proportional to the 
electric field, is much smaller than the characteristic relaxation rate in a crystal. The 
situation has changed in 1970, when Esaki and Tsu reported first measurements of the 
Ohm current in semiconductor superlattices [2|]. It was found that with increase of the 
electric field the current reaches some maximum value and then decreases - the phenomenon 
known nowadays as the negative differential conductivity. This pioneering work initiated 
experimental and theoretical studies of ordinary conductivity in the strong field regime, 
where the electric field cannot be considered as a perturbation. We note that, besides 
semiconductor superlattices, the negative differential conductivity was also observed for 
cold neutral atoms in (quasi) ID optical lattices subjected to a static force j^. A great 
advantage of the latter system over semiconductor superlattices is full experimental control 
over relaxation processes. Because of this control one can study with cold atoms both 
the Hamiltonian and dissipative dynamics of the carriers, i.e, Bloch oscillations and Ohmic 
current. 

A different yet related transport problem is the Hall current in quantum dot and antidot 
arrays (see, {4- 7|, to cite a few of hundreds relevant papers). These systems can be considered 
as 2D superlattices, where the 2D electron gas is subjected to a periodic potential with period 
of few hundreds nanometers. Here researches mainly focus on the effect of a magnetic field 
rather than on Bloch oscillations js]. Indeed, because of a large superlattice period already a 
modest magnetic field ~ IT causes dramatic changes in the electron energy spectrum. This 
new spectrum is then substituted in the Kubo-type equation to calculate the conductivity 
tensor aij = (Jij{B). During the last two decades many exciting theoretical and experimental 
results on the Hall conductivity in quantum-dot arrays were reported. These studies revealed 
very nontrivial dependence of the conductivity tensor on the magnetic field and the Fermi 
energy. Let us also note that in a typical quantum-dot array the Fermi energy of the 2D 
electron gas is in a Bloch (mini)band with large index n > 10, so that the single-band 
approximation, which is known to result in the famous Hofstadter's butterfly spectrum {q], 
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is not applicable. Because of this complication of the problem almost all theoretical studies 
on the Hall conductivity in quantum-dot arrays have been restricted by the linear response 
regime [lO|. 

Recently it has been noticed in Ref. [l3^] that the Hofstadter butterfly can be realized with 
cold atoms in a 2D optical lattice by introducing an artificial magnetic field, which imposes 
the Lorentz force on moving neutral atoms. To create this artificial field the authors of the 
cited paper suggested to use Raman transitions between the internal atomic states. This idea 
of an artificial magnetic field was developed further in other publications and today there 
are several different theoretical proposals for creating effective magnetic fields for atoms in 
a lattice (see, for example, the recent paper [IJ] and references therein) and a successful 



15| . Thus this is a matter 



realization of one of schemes for atoms in a harmonic trap 
of time when experimentalists will be able to mimic magnetic fields for cold atoms in a 
lattice, in the same manner as they now routinely mimic electric fields. Since optical lattices 
are much cleaner and controllable systems than quantum-dot arrays, this will open wide 
perspectives for studying the Hall conductivity in different parameter regimes, including 
the nonlinear response regime. The present work precedes these anticipated experimental 
studies. It presents a semi-analytical method for calculating the Hall current in the single- 
band approximation which, as mentioned above, can be easily justified for cold atoms in 
optical lattices. The method treats magnetic and electric fields on equal footing and, hence, 
goes beyond the linear response theory. In this sense we extend the Esaki-Tsu approach to 
the Ohm current in ID lattices onto the Hall current in 2D lattices. 



II. ESAKI-TSU DEPENDENCE FOR THE OHM CURRENT 

First we recall the reader few results on ordinary conductivity in the non-perturbative 
regime. To explain the negative differential conductactivity in semiconductor superlattices 
Esaki and Tsu used a kind of semiclassical approach, which resulted in the following depen- 
dence for the Ohm current: 

In Eq. ([1]) F is the electric field, d the lattice period, e the charge, 7 the relaxation constant, 
and the pre-factor /(T) accounts for the temperature dependence of the current [/(O) = 1]. 
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A microscopic derivation of the Esaki-Tsu dependence ([T]) was given by Minot in 2004 |16 |. 
This was obtained by solving the master equation for the carriers one-particle density matrix 

In this equation H is the single-particle Hamiltonian of a carrier in a biased superlattice, 

H = Ho + edFY,\l)l{l\, ]/o = -^^(|/ + l)(/| + /,.c.) (3) 
I ^ I 

and po is the equilibrium density matrix for F = 0, 

po ~ exp{-0o) , /3 = l/knT . (4) 

The model (EM]) results in Eq. ([T]) with the correct pre-factor f(T) = Xi(/3J)/Xo(/3J). 

The above microscopic derivation of the Esaki-Tsu dependence was revisited in Ref. 
with respect to the problem of ordinary conductivity with cold atoms in ID optical lattices. 
A weak point of the master equation ([2]) is that it is not in the Lindblad form. (Exclusions 
are the cases of zero and infinite temperature, where it can be rewritten in the Lindblad 
form.) Because of this drawback it may give wrong result for the velocity distribution of the 
carriers. However, it was confirmed that it gives qualitatively correct result for the mean 
velocity, i.e., the current. The master equation ([2]) will be our theoretical framework in 
studying the Hall conductivity in 2D lattices. 



III. THE MODEL 



We consider a quantum particle in a square lattice of side d in the x — y plane. The 
particle is subjected to an in-plane electric field F, aligned with the y axis, and a magnetic 
field B normal to the x — y plane. Using the tight-binding approximation and the gauge 
A = B{—y, 0, 0) for the vector potential the particle Hamiltonian reads 

H = Hp + edF \l, m)m{l, m\ , (5) 

l,m 

^0 = -y E (l^ + 1> ^) m|e^2™r. ^ _ :| ^ (|/^ ^ + i) (/^ ^| + . (6) 

l,Tn l,m 

The dimensionless parameter a in (|6]) is the Peierls phase, which is given by number of 
magnetic flux quanta per unit-cell area, a = eBd'^/hc. Besides the Peierls phase a and 
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the Bloch frequency ub = edF/h the other important characteristics of the system are 
the carrier effective mass, M* = cP{JxJyY^'^ /h^ the cyclotron frequency, ujc = eB/cM* = 
27ra{JxJyY^'^ /h, and the drift velocity v* = ceF/B = cPeF/ha. We note that for a charged 
particle (electron in a solid crystal) the Hamiltonian (jH]) is justified only in the limit of small 
a, where the cyclotron radius of the classical orbit essentially exceeds the lattice period. This 
is, however, not the case for cold atoms in optical lattices, where the actual parameter of the 
system is the Peierls phase but not the magnitude of a magnetic field [isl, Q- Hence, we 
impose no limitations on a and, without any loss of generality, one may consider |a| < 1/2. 
Our aim is to calculate the Hall {vx) and the Ohm {vy) currents, 

V^^y = Tllv^^ypst] , (7) 

where pst is the stationary solution of the master equation (|2]) and Vx,y the current operators, 

Vx = -^[Ho,x], X = d'^\l,m)l{l,m\ , (8) 

l,m 

and for Vy one has a similar expression. Substituting ([6]) in ([8]) we have 

*^ = ^E(l^ + l'^)(^He^''^""^-/^.c.) , vi,^^ = ^, (9) 

l,m 

and 

*. = ^E(l^^ + l)(^'H-/^-c.) , v^'^ = ^^. (10) 



IV. LANDAU-STARK STATES 

We shall perform calculations in the basis of the Landau-Stark states, which are the 
eigenstates of the Hamiltonian ([5]). To simplify equations, from now on we set the lattice 
period d and the Planck constant h to unity. 

One finds the Landau-Stark states semi-analytically by using the following ansatz jisl : 

^Z,m = -y^bm ■ (11) 

In Eq. ( ITT]) K = 2TTk/Lx is the quasimomentum, < n < 27r, is the lattice size in the 
X direction and we eventually let tend to infinity. Substituting (fTTj) into the stationary 
Schrodinger equation with the Hamiltonian ([5]), 



we reduce it to the following ID eigenvalue problem: 

- y + bm-i) + [Fm - J.^ cos(27ram - K)]bm = Ebm ■ (13) 

Equation f|T3l) is a kind of the ID Wannier-Stark problem and can be easily solved nu- 
merically. Labeling the solution by the discrete index u and scanning over the quasi- 
momentum n we find the energy spectrum E = Eu{k,) and the Landau-Stark states 

The energy spectrum and properties of the Landau-Stark states were studied in some 



detail in our recent work 19| devoted to the Hamiltonian dynamics of the system ([5]) . As an 
example, Fig. [1] shows the energy spectrum E = E^{k) for = Jy = 1, a = 1/10, and two 
different values of F. This figure is aimed to illustrate a qualitative change in the spectrum, 
which takes place around 

F,r = 27raJ, . (14) 

Namely, for F < Fcr the energy form a pattern with straight lines. The Landau-Stark states 
belonging to these lines are the transporting states, which transport the quantum particle 
in orthogonal to the field direction with the drift velocity f*, 

V* = F/27ia . (15) 



Having the Landau-Stark states obtained we calculate the current operators ([9]) and ( JTOjl 
in this basis. We have 

= '5(/t-/t0^xE^mH«^)^'L''n«^)sin(27ram-/€) , (16) 

m 

and 

{^.,.\%\^uw) = - '^O^ E [tU^^) - htU^)] C')(«:) = 5(«: - ^')v!^k^) . (17) 
Because of the presence the 5-function in ( IT61IT71) . Eq. (JTj) for the Hall current simplifies as 

1 /"^TT 



/ "d/^Tr[V^("•^)(/^)7^(^*)(/^)] , (18) 
Jo 

where 'RS^^\k) is the K-specific stationary density matrix. 



ZTT JO 



7^i5(/^) = {^.APst\'^u',.) = —TT^ril^-^rj-r^n^Z'i^) . (19) 

•J + z[E^'[k) - F^[k)\ 
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FIG. 1: A fragment of the energy spectrum of the Landau-Stark states for = Jy = oi = 1/10, 
and F = 1 (left) and F = 0.3 (right). The figure is borrowed from Ref. [l^ . 

V. LANDAU STATES 

Next we specify the equilibrium density matrix po- To have tractable results we shall 
consider the case where only the lowest Landau levels are populated. Thus we assume 

Po = nj7El$.)($.l, (20) 

where M = LyL^a and are the lowest energy Landau states. The density matrix (l20l) 
corresponds to Af fermionic carriers at zero temperature. Alternative, it may be considered 
as a density matrix of non-interacting bosons at a finite temperature. In what follows we 
adopt the latter point of view, where the relevant temperature interval is discussed in the last 
paragraph of the section. We would like to stress that our choice of the equilibrium density 
matrix is exclusively for the sake of easy interpretation of numerical results. In principle, 
one can consider an arbitrary po- This way the reported below results can be generalized to 
arbitrary temperature and arbitrary carrier statistics. 

Similar to the case of Landau-Stark states, one finds the Landau states semi-analytically 
by using the substitution (|TT|) . which reduces the eigenvalue problem for the Hamiltonian 
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FIG. 2: Integrated density of states of the system ([6]) for a = 1/20 (left) and a = 1/5 (right). The 
other parameters are Jx = Jy = ^ and Lx = Ly = 40. 

Hq to that for Harper's Hamiltonian, 

- yl&m+i + hm-i) - Jx cos{2nam - K)bm = Ebm ■ (21) 

Figure [2] shows the integrated density of states of the system ([6]) calculated on the basis of 
Eq. fl2T]) for a = 1/20 and a = 1/5. The states in fl20|) are associated with the first step 
in the integrated density. It is also easy to show that the length of this step is approximately 
given by the cyclotron energy huc = 2Tia{JxJyy^'^ ■ Thus our condition on the temperature 
reads fc^T ^ huc- At the same time, to have equal populations of the lowest Landau states, 
we assume ksT ^ A, where A is the width of the lowest magnetic band {2^. As seen in 
Fig. [2] the required condition A ^ ksT <^ TiUc is easier to satisfy for small a. 



VI. NUMERICAL PROCEDURE AND RESULTS 



The numerical procedure is as follows. We fix Lx and Ly and calculate the Landau and 
Landau-Stark states. The lattice size Lx defines the discrete step for the quasimomentum, 
which should be small enough to resolve main quasi-crossings in the energy spectrum in 
Fig. [H The lattice size Ly is arbitrary yet, to reduce the boundary effect when solving 
(fT3ll . Ly ^ 1/a. Next we calculate the K-specific matrices of the current operators and the 
stationary density matrix (fT9l) . We note that for a rational a = r/q the infinite matrix of 
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FIG. 3: The Hall (left) and Ohm (right) current as the function of electric field F for different 
relaxation rates 7. The other parameters are Jx = Jy = I and a = 1/10, the lattice size = 
Ly = 40. The straight dash-dotted lines are predictions of the linear response theory. Additional 
dashed line in the left panel shows the Hall current for 7 = 0, where the Ohm current vanishes. 

the current operators as well as K-specific density matrix obey the translational symmetry, 



(22) 



which further facilitates the numerical procedure. Finally, substituting these matrices into 
(fT8|) and integrating over the quasimomentum k we calculate the Hall and Ohm currents. 

The left panel in Figj3] shows the Hall current Vx as the function of the applied field F for 
a = 1/10 and different values of the relaxation constant 7. We begin with considering the 



case 7 = (dashed_ 
As shown in Ref. 



ine), which corresponds to the Hamiltonian dynamics of the carriers. 
19| . for the specified initial conditions (population of the ground Landau 
states) a weak static field transports the carriers in the orthogonal direction with the drift 
velocity f lT5|) . Thus in the weak field regime the dependence Vx = Vx{F) is approximated by 
Vx = F/2Txa. The increase in the Hall current continues till F reaches some critical value 
F*, where the function Vx = Vx{F) has the global maximum. An estimate for F* is provided 
by Eq. f ll4p . although we found F* to be systematically larger than F^r by a numerical factor 
1 < 2; < 2 [see Fig. Hl^a) below] . With further growth of the electric field we enter the regime 
of negative differential conductivity, where the Hall current decreases with increase of F. 

As mentioned above, one finds an explanation for the transition from positive to negative 
differential conductivity regimes in structural changes of the Landau-Stark states, which 
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FIG. 4: The Hall (left) and Ohm (right) current as the function of electric field F for 7 = 0.1 and 
different a. Additional dashed line in the right panel shows the Ohm current for a = 0, where the 
Hall current vanishes. 



take place around Fcr 19|]. We also give another explanation, which is based on the Landau 
states picture. Namely, using the Kramers-Hennenberger transformation an electric field is 
seen as periodic driving of the system with the Bloch frequency ojb = F. When ojb matches 
the energy gap between the ground and the first magnetic band [i.e., the length of the 
first step in Fig. [21 approximately given by the cyclotron frequency Uc], the driving induces 
transitions between the Landau levels and we observe the local minimum in the dependence 

= Vx{F) which always precedes the global maximum. Thus the necessary condition for 
the negative differential conductivity regime can be also formulated as a requirement that 
the Bloch frequency exceeds the cyclotron frequency. 

The other (solid) lines in Fig. [3] show the dissipative Hall current for 7 = 0.1,0.5, 1. It 
is seen that a finite relaxation rate suppresses the Hall current and smoothes fine features 
of the dependence Vx = Vx{F) for the non-dissipative Hall current. In addition to Fig. [3t^a) 
figure m^a) shows the Hall current for the fixed 7 = 0.1 and different a. The vertical dashed 
lines in this figure indicate the critical electric field f|T^ for each case. 

The right panels in Fig. [3] and Fig. H] show the Ohm current. It is seen in Fig. [3] that 
larger relaxation rates suppress the Hall current but enhance the Ohm current. In the limit 
7—7-00 the Hall current vanishes and the dependence Vy = Vy{F) for the Ohm current 
approaches the Esaki-Tsu dependence ([I]). Alternatively, one recovers the Esaki-Tsu result 
by considering the limit a — )■ 0, see Fig. Hl^b). In Fig. [3] and Fig. H] we also depict predictions 
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of the linear response theory, v = crF, where the off-diagonal and diagonal elements of the 
conductivity tensor are given by the Drude-type formulas, 

and we approximate the cyclotron frequency by uJc = 2T{a{JxJyY^'^ ■ 



VII. CONCLUSIONS 

We considered the quantum particle in a 2D lattice subjected to (real or artificial) electric 
and magnetic fields and calculated the Hall and Ohm currents as functions of the electric field 
magnitude. Although the obtained dependence = Vx{F) for the Hall current resembles 
the Esaki-Tsu dependence Vy = Vy{F) for the Ohm current in the absence of magnetic field, 
the physics behind these two dependences is completely different. Indeed, the Esaki-Tsu 
dependence for the Ohm current appears due to an interplay between Bloch oscillations and 
relaxation processes and the Ohm current vanishes if 7 = 0. Conversely, the Hall current in 
the transverse direction takes place even in the absence of dissipation. The actual reason for 
the Esaki-Tsu like dependence for the Hall current is a qualitative change in the structure 
of the Landau-Stark states, which happens around Fcr- Note that the condition F = Fcr 
means the Bloch frequency to coincide with the cyclotron frequency. Thus the Ezaki-Tsu 
like dependence for the Hall current is the result of an interplay between Bloch and cyclotron 
oscillations but not Bloch oscillations and relaxation processes. 

Concluding the discussion we would like to stress that in this work we do not addresss the 
quantum Hall effect. The latter phenomenon occurs for fermionic carriers when the magnetic 
field or the Fermi energy are varied. It would be interesting to study the quantum Hall effect 
in the non-perturbative regime, where the conductivity tensor depends on the electric field 
magnitude. This problem has been addressed in already cited papers [lliQ. However, these 
papers analyze the Hall conductivity specifically with respect to semiconductor structures 
and employ a different model. 
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